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Abstract 



>< 
;_( ■ The one loop effective action in quantum field theory can be expressed as 

a quantum mechanical path integral over world lines, with internal symme- 
tries represented by Grassmanian variables. In this paper, we develop a real 
time, many body, world line formalism for the one loop effective action. In 
particular, we study hot QCD and obtain the classical transport equations 
which, as Litim and Manuel have shown, reduce in the appropriate limit to 
the non-Abelian Boltzmann-Langevin equation first obtained by Bodeker. 
In the Vlasov limit, the classical kinetic equations are those that correspond 
to the hard thermal loop effective action. We also discuss the imaginary time 
world line formalism for a hot 0^ theory, and elucidate its relation to classical 
transport theory. 



1 Introduction 

In classical kinetic theory, a covariant formalism can be obtained in terms 
of phase space averages over the trajectories of particle world lines ^. In 
quantum field theory, there are many instances at finite temperature and 
density where classical ideas are relevant, and where a classical kinetic picture 
would be useful. However, it is not immediately apparent how one recovers 
the classical world line picture directly from quantum field theory. This is 
especially problematic in theories with internal symmetries. 

Fortunately, in the last decade, there has been a considerable body of 
work relating the one loop effective action in quantum field theory to quan- 
tum mechanical path integrals over point particle LagrangiansQ. For a survey 
of recent developments, see the review by Schubert @]. These recent devel- 
opments follow from the key insight by Berezin and Marinov 0] that internal 
symmetries such as color and spin had classical analogues in terms of Grass- 
manian variables. These obeyed classical commutation relations which, when 
quantized, gave the usual commutation relations for spin and color. 

Brink, DiVecchia, and Howe used these Grassmanian variables to con- 
struct a classical Lagrangian for spinning world lines in an Abelian back- 
ground field p. Subsequently, Balachandran et al. and Barducci et al. [Q 
wrote down the following Lagrangian for a classical colored particle in a 
non- Abelian background fieldQ: 



L = -mJxf.x'^ + i\\Dabh ■ (1) 



Here the Aa(r) with a = 1,---,N are Grassmanian dynamical variables, 
and DabK = Aa + igx^A^T^f^Xf,. The variables T^'s are N x N matrices in 
an irreducible representation of the Lie algebra of the group. The Euler- 
Lagrange equations of motion are deduced in the usual way from the above 
Lagrangian. It was shown by Balachandran et al. and by Barducci et al. that 
these equations are precisely the equations written down nearly thirty years 
agoby S.K. Wong §. 

The connection of the work on point particle Lagrangians to quantum 
field theory was first made by Strassler [^. He showed that the one loop 
effective action in quantum field theory could be expressed in terms of a 
quantum mechanical path integral over a point particle Lagrangian. For an 
Abelian gauge theory, Strassler showed that this Lagrangian was identical 



^The connection between fields and particles is of course relatively ancient. It goes back 
to the works of Feynman and Schwinger |g]. Also, Polyakov in his book g demonstrates 
the relation between particle Green functions and the n-point Feynman amplitudes of 
quantum field theory. 

^They also considered the case of a classical colored, spinning, particle in a non- Abelian 
background field. We will not discuss this case here. 



to the one written down previously by Brink, DiVecchia, and Howe [|1^ . 
For non-Abelian gauge theories, Strassler wrote down an expression, but it 
did not have a Lagrangian interpretation since the trace over color-the Wil- 
son loop-was kept explicitly. The corresponding world line Lagrangian was 
first obtained by D'Hoker and Gagne using a coherent state formalism [p 



Remarkably, as observed by Pisarski and Tytgat [T^, the point particle La- 
grangian D'Hoker and Gagne obtained by integrating out fermions coupled 
to a vector background field, coincides with the Lagrangian in Eq. (|^). 

Pisarski and Tytgat have noted that since one can write the one loop 
effective action as a path integral over world lines, it is a possible explanation 
of the apparently mysterious result that classical kinetic theory in terms of 
the single particle distributions f{x,p,Q) gives rise to the well known HTL 



effective action for finite temperature gauge theories |jT2[. In this work, we 
extend the previous work on the vacuum world line formalism, and develop 
a real time many body formalism, which may be applied to a wide range 
of many body problems at finite temperature and density. In particular, we 
will focus our attention here on the apparently mysterious results that were 



discussed in Ref. |jT2|, and show how they may be understood in the many 
body world line formalism. Before we can be more explicit, we need therefore 
to discuss the status of recent work in finite temperature gauge theories. 

At very high temperatures, the physics of soft modes in gauge theories 
may be described by effective field theories ||13| . Integrating out modes with 



momenta p ~ T, one obtains an effective theory at the scale p ~ gT, charac- 
terized by a Debye screening mass m|, oc g^T"^. The gauge invariant effective 
action describing these modes is often referred to as the Hard Thermal Loop 



(HTL) effective action |1I4| , |T5|, |T6[. Recently, there has been much progress in 
describing very soft magnetic modes at the scale p ~ g^T by systematically 



integrating out modes at the scales T and gT [|r^, |T^, |T^, |2y, |21[] 



It was first shown by Blaizot and lancu [^| that the non-local hard 



thermal loop effective action could be derived in a local kinetic approach 
by systematically truncating the Schwinger- Dyson equations. The kinetic 
approach was very useful since it lead to an intuitive picture in which soft 
modes p ~ gT satisfy classical equations of motion and are coupled to hard 
modes p r^ T . These hard modes in turn satisfy the collisionless Boltzmann 
equation. 

This Schwinger- Dyson approach is a perfectly reasonable way to proceed. 
However, since the approximations are made at the level of the operator 
equations of motion, it can become cumbersome. A local effective action was 
developed by lancu |2^, which is more useful, and indeed has been the basis 



of some recent work [0, |J] . 

Some time later, Kelly et al. P3| wrote down an alternative transport the- 



ory for classical colored particles. These particles obey the classical equations 
of motion of a spinless colored particle, coupled to a non-Abelian background 
field y4°, that were first written down by Wong 0. 

In the picture of Kelly et al., the single particle distributions f{x,p,Q) 
are defined for an extended phase space of classical colored charges. (In this 
regard, their approach is similar to the earlier work of Heinz p6[ and col- 



laborators on a classical transport theory for colored charges.) They showed 
that this distribution obeys the collisionless Boltzmann equation. Making 
suitable approximations to their kinetic equations, they showed that they 
also recovered the HTL effective action. 

The classical transport theory of Kelly et al. is very simple and elegant. 
However, the formalism is ad hoc, and it was not clear whether such a for- 
malism could be "derived" in some systematic approximation. It was also 
not clear whether the distributions / were distributions in the statistical 
mechanics sense, since no phase space averaging had been performed in the 
derivation. Finally, as pointed out by Blaizot and lancu [^ , it was not clear 
why color could be treated as a classical degree of freedom. 

The problematic issue of statistical averaging in the work of Kelly et 
al. was clarified recently in the work of Litim and Manuel [1^. In this work. 



we show how, in the real time many body world line formalism, the classical 
transport theory of Kelly et al., as improved by Litim and Manuel, may be 
obtained from the one loop effective action in QCD. We obtain the same 
set of transport equations they do. Since our results are derived from the 
QCD one loop effective action, one can, in principle, go further. We will not 
attempt to do so in this paper, but save that effort for later work. 

This paper is organized as follows. In section 2, we summarize the world 
line formalism for the one loop effective action in QCD at zero temperature. 
In section 3, we discuss the point particle Lagrangian that results from the 
one loop effective action. We show that Wong's equations follow from this 
Lagrangian. We comment on the classical commutation relation among the 
various dynamical variables and their quantum counterparts. Next, in section 
4, a many body world line formalism is developed, and applied to the one loop 
effective action. In the following section, it is shown that the saddle point of 
this effective action gives the finite temperature classical transport theory of 
Litim and Manuel. As shown by these authors, their results agree with those 
obtained previously by Bodeker |l^ for soft modes in hot gauge theories. The 



final section explores the implications of our result, and discusses further 
avenues of research. In appendix A, we work out, in the imaginary time 
approach, a simple example illustrating the finite temperature world line 
formalism. We also point out, in this context, the connections to classical 
transport theory. 



2 The one loop effective action and the world 
line formalism in quantum field theory 

In recent years, the world line formalism has become popular as a technique to 
compute one loop scattering amplitudes in gauge theories. Initial interest in 
such techniques followed from the work of Bern and Kosower |]28|, who derived 
rules from string theory for computing one loop scattering amplitudes in 
gauge theories. Strassler showed that similar rules could be obtained directly 
in gauge theories by writing the one loop effective action, in background 
field gauge, as a one dimensional quantum mechanical path integral over the 



"world line" of a point particle in an external background field |10 



We consider, with Strassler as a guide, the example of the Lagrangian of 
a massless scalar field coupled to a background Abelian gauge field A^{x), 



where D^ = d^ — igA^. The one loop effective action is 



T[A] = -\og\det{-D^ 



Tr \ok(-D' 



(2) 



(3) 



(The general case of Dirac spinors coupled to scalars, vectors, axial scalars 
and vectors, and tensors, has been considered by several authors^) 
We then use the trick 



log(a) 



Ji y Ji Jo Jo t ^ 



dy 
'I y 

to re-write Eq. (^ as 

r°o dT f d^p 



°o dt 



(4) 



T[A] 



T 



(2vr; 



4 {p\ exp ( -^eT ij) + gA{x)f ) \p) 



/ -—Af VxVexp - dr l—x'^ + igA[x{T)]-x 
Jo T J Jo \2e 



(5) 



where V denotes the path-ordering. 

Above, e is the einbein (the square root of the one dimensional metric- 
taken here to be an arbitrary constant) and 



M = I Vp exp I — - / dr ep'^ 



for a compact review, see the talk by C. Schubert |2£ 



is a normalization constant. Analytically continuing this expression to Minkowski 
space-time and letting e — > —e, one obtains in the general case of A^ being 
a matrix in the group representation R of the scalar field 

r°° dT r [ /"^ / 1 

T[A] = / —^fVx^TRVexpi dr i—x'^ - igA[x{T)] ■ x 

where Tr^^ denotes the trace over the internal states. This expression can 
be factorized into a product of exponentials-the A dependent piece being a 
Wilson loop. It is therefore the expectation value of a Wilson loop of the 
background field, in a particular ensemble of loops. 

Strassler showed that one particle irreducible (IPI) Feynman diagrams, 
with A^ external background gluons and a scalar loop, could be constructed 
by expanding the above expression for the one loop effective action to order 
g'^ . Expanding the gauge field in a set of plane waves with definite mo- 
mentum, polarization, and color charge, he re-wrote the above IPI-action as 
factorized products of the traces over color and the N proper time integrals 
(and permutations thereoff). He then considered several examples-and com- 
puted, for instance, the gluon vacuum polarization in a theory with Dirac 
fermions and complex adjoint scalars. 

Our objective here is slightly different. We wish to write Eq. (^ in 
terms of a point particle Lagrangian of scalar particles which carry internal 
symmetries and are coupled to an external field. (The extension to this work 
to fermions is straightforward but will not be discussed here.) Towards this 
end, we shall here follow the later work of D'Hoker and Gagne (see also 



the related work of Mondragon et al. [0). D'Hoker and Gagne proved the 



following identity for an n x n Hermitean matrix M{t), 

Tr V exp (t J^ dr Mir)) = (^) " E /^^ ^^^^^ g^'/'CAt A+f-i) 

X exp ( - r dr[AU - ^A^MA] j . (7) 

Here = 2TTk/n, where k = 1,- ■ ■ ,n. As noted in the introduction, the 
A^ and A are independent dynamical Grassmanian variables (the subscript 
"AP" denotes anti-periodic boundary conditions) which act as eigenvalues of 
creation and annihilation Fermi operators, respectively, generating arbitrary 
finite dimensional representations of the SU{Nc) symmetry group. These 
are discussed further in the following section. Also, the value of A^A can be 
evaluated at any arbitrary value of r. 

In the adjoint case, the matrix M stands for the adjoint gauge field A'^T'^, 
with n = Nl — 1. Substituting the above identity in Eq. (H), we can write 
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the one loop effective action for the Lagrangian in Eq. (^ as, 

^27rJfc(AtA+f-l)/n 



(8) 



2 
fc = l 




with 



. 2 

^p{^) = ^ + AU + igx^X'^A^X . (9) 



The important point to note here [T^ is that the above Lagrangian, contained 



in the one loop effective action, is precisely the point particle Lagrangian 
of Balachandran et al., and of Barducci et al. in Eq. (0). This result is 
extremely suggestive, and is indeed the starting point for our work. The 
above result in Eqs. (H) and (|^) can be easily extended to the full QCD case- 
the expression for T[A] above will also include a path integral over classical 
spins (see Ref. ||11|). Again, this Lagrangian is identical to the generalized 
Wong Lagrangian 0. 

In the following section, we will discuss Wong's equations and show how 
they follow from Eq. (H). 

3 Wong's equations and the world line La- 
grangian 

In recent years, Wong's equations for classical charged particles interacting 
with classical non-Abelian gauge fields have received a considerable amount 
of attention in finite temperature applications |]25|, |2^, ^. Wong wrote down 
the following set of equations 

m— -^ = mvj = p^ (10) 



gviTi{QjF^^{xj)) (11) 



dr 

(¥[ 

dr 

Qi = -z(?[q,,<A'^] (12) 

D'F,^ = 3^ (13) 



and 



JMia^j 



gY. fdTQr{T)vl{T)6'[x-xj{T)] . (14) 

1=1-^ 



Above, r is the proper time, Q\t) = QHt) Ta is the color charge of the J-th 
particle in the adjoint representation, and {D^Y'^ = 6"'^df^ + gf°'^'^A'i is the 



usual covariant derivative. The second equation is the generahzation of the 
Lorentz force equation, the third equation describes color precession. The 
fourth equation is the Yang-Mills equation, which gives us 

D,f = , (15) 

where j^ is the net color current generated by the classical system of K 
particles. A similar set of equations including the spin was also obtained by 
Heinz |3^ starting from the Dirac equation for quarks. 

Balachandran et al. showed that the above set of equations can be ob- 
tained from the point particle Lagrangian in Eq. (P). We will summarize 
their results below. The action 

S = -jd'x ^F^^F^'^'^ + j dr C,{r) , (16) 

is a) real upto a total time derivative, b) re-parametrization invariant, under 
r — » r', and c) invariant under gauge transformations in which both A° and 
Aq are transformed. Thus under an infinitesimal gauge transformation ^°, 

^M — ^ {AiY = Ai-gU,ePAi-d,e\ 

Xa — > {\a)' = Xa + ^g9''T:,9b. (17) 

Redefining A^ — > iX^ , one obtains the equations of motion 

Xa + tgx,A':T:,X, = 0, (18) 



and 



where 



d ( X 



' ' + gQ^'F^x = , (19) 



Q° = XiT^^X, . (20) 



From the preceding equation, and from Eq. (p!8D , we obtain the equation of 
motion for Q" in Eq. (0). 

Finally, taking the functional derivative of Eq. (p!6D with respect to yl^(x), 
we obtain the Yang- Mills equations in Eq. ([13D , with the current j^ defined 
as in Eq. (^. 

In the Hamiltonian formalism, x'^, Xa and X^ are independent dynamical 
variables. Their momentum conjugates are 

P, = '^-gQ'^Al. 

Pa = ^Xi. 

Pi = 0. (21) 



respectively. Due to the second class constraints, we have to introduce Dirac 
brackets. For two dynamical variables A and B which are even elements of 
the Grassman algebra we have, 

{^,B}„ = ^f^ a 5 9 _^5 9 _^S aU, ,22, 

^ • ID yg.j:^gpf dPi-dxi' dXldK dKdXi) ^ ' 

This gives rise to the following Dirac brackets, 

{Q'^^Q^}^ = U,Q'- (23) 

Upon quantization, these give the usual Heisenberg commutation relation 
for P and x, and |Aa,AU = 5ab- Also, \Q",Q^\ goes over into the usual 
commutator. As pointed out by Balachandran et al. and by Barducci et al., 
the quantized A^'s and A's act as Fermi creation and annihilation operators, 
and generate various representations of the color group. The n creation and 
annihilation operators span the Clifford algebra C2n- The sum over (/)'s in 
Eq. (^ may therefore be understood as restricting the states generated by A 
and A^ to a particular irreducible representation of the group (see for instance 
Eq. (3.9) in Ref. §). 

We have dwelt at some length on the classical and quantum commutation 
relations among the dynamical variables because they will be useful later on 
in understanding the transition from the one loop effective action to classical 
transport theory. 

Before we end this section, we should mention that the generalization of 
the action in Eq. (|1^) to include classical spins was also performed by Bal- 
achandran et al., and by Barducci et al. The classical spins are also Grassma- 
nian variables, and upon quantization, go over into the Dirac 7-matrices. One 
obtains the generalized Wong equations, including a generalized Bargmann- 
Michel-Telegdi equation for the time evolution of the Pauli-Lubanski spin 
four-vector. The generalized point particle Lagrangian is again identical to 
the one that appears in the one loop effective action for fermions in a classical 
background field. 

4 The many body world line formalism 

In this section, we will extend the formalism of section 2 to study a many 
body system of world line scalars. This fills in the step where one goes from a 
quantum field theoretic description in terms of ensemble averages over fields, 

9 



to the classical transport approach where one takes phase space averages 
over world line trajectories. These averages will of course correspond to the 
multi-particle distributions of classical transport theory. 

In general, for the path integral of a many particle system, one has an en- 
semble of initial and final conditions0( which do not necessarily coincide!). In 
addition, not all initial and final conditions are weighted equally-the prob- 
ability for these being given by a density matrix. For example, at finite 
temperature, in the real time formalism, the initial conditions are specified 
at a fixed initial time in the complex t-plane, the probability being given by a 
Boltzmann weight. This is reasonable because in a thermal system operators 
are of course not "sandwiched" between vacuum states. 

We begin again with the full path integral for a massless scalar coupled 
to an external classical gauge field A^. The generating functional is 

itjj] exp (i^^^D^^jjj 

= j[dilJi] [#/] /C[^,] Z[A^] Z*[A2] (24) 

where the subscript C represents the usual Schwinger-Keldysh closed time 
path. One may consider the two pieces of real time paths as representing the 
two time evolution operators needed to evolve a density operator. Here, Ai 
is the background field in the positive time direction, A2 is the field in the 
negative time direction and 



Z[A] = / [#i] exp (^V^JDVi) (25) 



The density operator itself is represented here by the matrix element 

/C[7/^,] = (V^ilAnitlV^^) (26) 

where ipi are eigenvectors diagonalizing the Hermitan matrix Pinit- 

To rewrite the above in terms of determinants, we note that Z[A\ can be 
written as 

Z[A]-^ = n ^n (27) 

n 

where 6'„ is the eigenvalue of the equation 

D%,, = eniJn (28) 

^For an excellent reference on path integrals for many particle systems, we recommend 
the book by Negele and Orland p3| . 
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with the boundary conditions given by V^(tinit) = i'i and '^/'(tfin) = "ipf- 
Taking the logarithm of Z[A] yields 

-\ogZ[A] = Y.^oge^ = J2{n\logD^\n) 

n n 

= Y. f d'^xd'^y J2{n\y,\') {y,\'\\ogD'^\x,\) {x,\\n) 

A,A' 

= J2 [ d'xd'yJ2^Pliy,X'){y,X'\\ogD'\x,~X)M^,>^) 

A,A' 

(29) 

Here A is the color label of a particular representation of the algebra^. In 
contrast to the vacuum case, one cannot simply say here that 

J2£{y,x')Mx,~^) = s^'Hx-y)h'x 

n 

since z/'n can only span functions with the same boundary conditions. In the 
vacuum case, one is able to do this because boundary conditions at t = ±cx) 
are usually trivially the same for any reasonable observables. We now use 
the identity in Eq. (^, and re- write Eq. ( [25| ) using the same line of reasoning 
as in section 2. (For details of this procedure, see Eqs. (3.16)-(3.22) in the 
first paper of Ref. |TT|). We obtain 

Y\A,i\ = -log Z[A,^] 

d'xd'yJ2ax,y,X,X') -^^ - ^W ^AtPA 
f^, JO I Jx,x Jap 

A, A 

X {j^yit e'-'^'''^^-'^/-exp (^^l'^drC,{z,X,X',A)^ ) . 

(30) 
We have replaced in the preceding equation 

^ i;i{y, X') Mx, A) = e(x, y, A, A') , (31) 

n 

where ^ is a functional of both ipi and ipf- Also, Cp is identical to the point 
particle Lagrangian in Eq. (|^). 



^As in the earlier sections, the A's symbohze the eigenvalues of the coherent state 
creation and annihilation operators which are generators of arbitrary finite dimensional 
representations of internal symmetries. We can write ipniv,^^) = i'nix) + Xl^tpn,a{x) + 
>^l>^liJ7i,ab{x) • • • as in Ref. 0. 
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Using Eq. (p^ ) and Eq. (pOD above, one obtains the following many body 
path integral, 

Z = j[di\ exp (-Gie]) ^ [dA] exp (^S-eff) , (32) 

where 

5eff [A, e] = - / ^'a; \f;,F^^^'' + r[A, e] . (33) 

Here we have made the assumption that the unknown K, is such that it 
permits Eq. ([3^ ) with a positive definite G[^]. We do know from Eq. (^) 
that /C is real since pinit is Hermitean. What we do not know is whether the 
Jacobian associated with the variable change (^ipi,ipf) ^ ^ is real. This is 
not clear since it is not guaranteed that ^(x, y, A, A') is real. However we will 
show shortly that when the saddle point approximation is made to obtain 
classical world lines, only the real part contributes. And in that case, it is 
likely that the relevant part of G[Q is positive definite. 

The many body path integral in Eq. (|32|) thus includes an average over all 
possible ^ with the weight e"'^'^^ . Since the ^'s here are determined exclusively 
by solving the determinant for the hard modes one arrives at the following 
physical picture. The above effective action can be thought to represent the 
dynamics of soft modes which are coupled to a background of hard modes 
represented by the point particle world lines. That coupling is weighted by a 
functional G[^]. The above form of the effective action is reminiscent of the 



small X effective action |34|. There, by making a separation of scales between 



small X and large x, it was shown that G[^] obeyed a non-linear Wilsonian 



renormalization group equation [^. In finite temperature QCD, a similar 
separation of scales occurs at momenta p ~ T,gT,g'^T. Thus G[C,] ~ Ga[^], 
where A is the scale {gT <^ A <^T) separating hard and soft modes fl^, p7|] . 
It will be interesting to see if renormalization group techniques^, analogous 
to those employed in the small x effective action, can be used to determine 
the evolution of G[^] as a function of the cut-off at finite T |^ ^ . 

We should comment on the real time contour C in the path integral over 
the effective action in Eq. (0). In the real time formalism of quantum field 
theory, for thermal initial conditions, the path consists of a piece where the 
field theory is evolved in imaginary time from t = to t = —i(3 and two time 
slices from — cx) to -|-cx) and back from oo — if3 to — oo — i(3. Propagators 
have a 2 X 2 matrix structure, with the off-diagonal terms corresponding to a 
flux of on-shell, thermal particles. The corresponding problem of many body 



propagators in the world line formalism has been addressed by Mathur [^ 

^For a discussion of the real time renormalization group, see the book by Golden- 
feld Hi. 
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He considers a more general initial state than a thermal one, with the only 
restriction on the initial density matrix being the constraint that it can be 
represented as p = e^'*'^', where B is quadratic in the field operator 0. The 
thermal distribution is a particular instance of such 'exponential of quadratic' 
density matrices, which have the property that correlators computed with 
these matrices, satisfy Wick's theorem. 

In the first quantized language of world lines, one obtains the following 



picture [^. Since the world line particle moves only in one direction in 
proper time, the paths backwards in time correspond to switching the sign 
of the einbein discussed in section 2. Specifying the amplitude ^ to a re- 
verse orientation of proper time is equivalent to specifying an 'exponential of 
quadratic' density matrix. Replacing ^ by a ^-function, as discussed previ- 
ously, would let us recover the Feynman propagator. For the particular case 
of classical transport theory in hot QCD to be discussed in the next section, 
we will not need to specify the real time contour, but it will be relevant for 
future detailed computations. 

To summarize, in this section, a function ^ is obtained, and a path inte- 
gral introduced over all ^ with an unspecified weight G[^]. Effectively, what 
one has done is to introduce a density matrix p corresponding to the hard 
modes in the system. The density matrix of course is required to satisfy the 
condition 

Tr(p) = /Me]exp(-G[e]) = l. (34) 

5 From the QCD one loop effective action to 
classical transport theory 

In the previous section, we extended the world line formalism discussed in 
section 2 to a many body context. Our objective in this section is to derive 
from the one loop effective action in many body QCD under suitable ap- 
proximations, the classical transport theory of soft classical modes coupled 
to the world lines of the hard modes. In particular, we will apply the world 
line formalism to gauge theories at the very high temperatures where weak 
coupling techniques are applicable. We will see that we recover the coupled 
set of transport equations that were written down by Litim and Manuel [|l^ . 
It was shown by these authors that their results, in turn, agreed with the 
earlier work, in an apparently very different approach by Bodeker [0. 

In hot QCD, the leading high T contribution to the one loop effective 
action for gluons is the same as the action for the example we discussed 
in previous sections-namely, that of a complex scalar field in the adjoint 
representation. We see this as follows [^. Working in the background field 
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gauge, and integrating out the fluctuations to quadratic order in the hard 
fields, one obtains 



Scs = Ssoft + Tr (log(-D2)) - ^Tr (log [-D^g^, + 2F^, 



(35) 



where D^ is the usual covariant derivative for the background field. The 
above equation will also have in general a gauge fixing term-with a gauge 
fixing parameter a. For simplicity, we have eliminated this term by setting 
Feynman gauge, a = 1. In Ref. [|l[], it was shown that the a dependence is 
suppressed by powers of 1/T. 

One can expand the above equation, writing 



5', 



eflf 



S, 



soft 



Tr(log(-D2)A _Tr log 



2F, 



j^lU 



D2 



S,,,, + T[A]. 



(36) 



By power counting arguments |T1], the terms containing F^'^ are suppressed 
by powers of T at high temperature. Hence, T[A\ for hot QCD is the same 
as the r we derived in the previous section. 

We now wish to derive the classical transport theory that follows from the 
above effective action. If we insert a functional derivative 5/5A in Eq. (|3^) , 
we have the identity 

j [di] exp (-G[e]) j^ [dA] -^ exp {iS,^) = . (37) 

This is of course equivalent to the operator equations of motion 

{D^F^'^ - r) = Tr (p [D^F^"" - J")) = . (38) 

From Eq. ( pT]) and Eq. (pOD, one obtains 



(^'"(x)) 



-cm 



J^[dA]e''^« JdrdyJ2U{r,y,~X,~X') 

A, A 



/o T ir,A Jap \tJ n ^ 



k=\ 



/o M'^'^fx) 



z — \ : — ^ o '^0 



(39) 
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Let's now take the saddle point of the Quantum Mechanical path integral 
in the expectation value of the currentQ At present, this step can only be 
justified a posteriori. However, it can be anticipated. One expects the world 
lines of the hard p ^ T modes to obey classical trajectories a la Wong since 
their coupling to the soft modes is weak at very high temperatures. (This 
also happens to be the case for the coupling between small x and large x 
modes in the small x effective action.) This information should of course be 
obtained self-consistently from the effective action. We hope to address the 
validity of this classical approximation at a later date. 

Taking the saddle point of the Quantum Mechanical path integral, we 
obtain 

(J^'"(x)) = J m exp {-G[Q J^ [dA] exp (zS^^) r^'J{x) , (40) 



where 



J^A^^) = Jdrdyj:ar,yXX')Af' f^[^y - 



X,X' 

X expii dr'Cpizc, Ac, Aj, A) 









Jo 


drg 


'XlT'^X^z", 


5^^\x 


^c 


■ir)) 


(41) 




In 


Eq.(40) 


the effective action is now 


defined to be 5'eff 


= s 


soft + Tel 


where 




Fci 


-f 


drdy > ]C{r, 

X,X' 


,y,X,X')Af' 


Jo T \t) 


/ 

exp 'i 

V 


Jo 


dr'Cp^z, 


:;) Ac, Ac, 


A) 






















(42) 



The subscripts c on the variables correspond to the classical paths which 
have the endpoints (r. A) and {y, A'). Note that for the classical path the r 
integral is really independent of T because rescaling r -^ Tt does not change 
anything. (Recall from section 3 that the classical point particle Lagrangian 
is reparametrization invariant.) Hence, extending the range of integration, 
< r < T — > < T < oo should not modify the result. Performing this 

^The sum X]fe=i e"^'^"^ A+n/2-i) j^^ ^j-^^^ path integral in Eq. (|39| ) is analogous to the 
well known GSO projection in string theory. As discussed in section 3 (after Eq. (P2|)), this 
projection restricts the world line paths to lie in a particular irreducible representation 
of the group. We will assume implicitly hereafter that this restriction is obeyed by the 
classical path. 
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operation results in 



J{,A \^) 



g jdrdyY. U{r, y, A, A') j^ dr Iq'^{t, r, y, A, A') i,^(r, r, y, A, A') 

X 6^'\x-z{T,r,y,~X,~X'))V\. (43) 

We have here used Eq. (^) to introduce the color charge Q°'. Also, we have 
redefined ^ to absorb the infinite normalization constant coming from A/"' and 
the T integral. A discussion of the renormalization of infinities in the first 
quantized formalism is outside the scope of this paper. We refer interested 



readers to the discussion in Ref. ||42|| . 



Now introduce the dummy integrals 

I dQ5^'^\Q-Q{T)) = l and f d''p5^''\p - p{r)) = 1 . (44) 

The color charge measure imposes the constraint (see footnote 7) that the 
Casimir invariants for the group are conserved. For SU(3) [p5|] , 

dQ = S^'>Q 6{QaQ^ - ga) S{dabcQ''Q'Q' - gs) , (45) 

where g2 and q^ fix the values of the Casimir invariants, and dabc are the 
symmetric constants of the group. The color charges spanning the phase 
space are dependent variables but, as argued by Kelly et al., they can be 
formally related to a set of independent phase-space Darboux variables [^ . 
We will not further address the issue of the color measure in this paper. 

Substituting the dummy integrals for p and Q, making use of Wong's 
equations, and switching the order of integration, the current can now be 
re-written as 

J^Ai^)=9 J d^pdQp>^Q^f^Ax,p,Q). (46) 

We have introduced a function f^^A{x,p, Q) which is defined by the relation 

f^^A{x,p,Q) = / dr drdy^U{r,y,\,\')6^'^\x- z{T,r,y,\,\')) 

° A,A' V 

X 5^-\Q-Q{T,r,y,~X,~X'))6^'\p-p{T,r,y,~X,~X'))y (47) 

Consider now what f^^A represents. First note that written in this way, it 
is clear that only the real part of ^ contributes to the color current. This can 

16 



be easily shown by noticing that ^(r, y, A, A') = ^*{y,r, X',X) and renaming 
r ^-> y, A <->• A'. Classical paths z and Q are not affected by this exchange 
because given the boundary conditions, these are fixed. The meaning of 
^(r, y, A, A') is that of the probability density of a particular path that begins 
with (r, A) and ends with {y, A'). However, since we are dealing with classical 
paths here, one can equally well specify a trajectory by the position and the 
velocity and the value of the charge at any fixed time|. (Note that since 
Eq. (|12]) is a first order differential equation, one needs only one boundary 
condition.) In particular, we can do it at r = 0. If we perform the trace 
explicitly, there will be Jacobian corresponding to the affine transformation. 
For a discussion of this point, we refer the reader to a paper by Brandt, 
Frenkel, and Taylor [Q-in particular Eq. (2.14) therein. 

Returning to Eq. (0), one can show that f^^A satisfies the equation]^ 



r^l^ 



dxf^ 



2gTriQF^n^ + ^9M[A„Q]^)] fU^,V,Q)=^- (48) 



Although it appears as such, the above equation is not quite the collisionless 
Vlasov equation, since f^^A here is a microscopic quantity and not a single 
particle distribution in the statistical mechanics sense. 

To obtain an equation for the single particle distributions, averages both 
over S, and A must be performed. Keeping in mind that the expressions here 
are functionals of A and .^, we have 

= p'^(^(4A(x,p,Q))-2^Tr(^(QFn^4A(x,p,g)^ 

- 2gTi U[Q,A^]^)fUx,p,Q)\\ (49) 

Here the bracket (■ ■ ■) indicates the functional averages over A and ^ with 
the weights specified in Eq. (^01). The corresponding equation for A is 



{D^F^ = {Jn (50) 



where (J) is expressed in terms of f^^A via Eqns. (1^) and (^q). Note that 
the left hand side of this equation consists of 2 and 3 point functions of A. 

The equations (|4^) and ([5^) are our main results in their most general 
form. They are not the Vlasov equations as noted by Litim and Manuel [^]. 
This is simply because these equations know about 2 and 3 point correlations, 
whereas in the Vlasov equation only the average values appear. Ignoring 
these correlations, one would get the Vlasov equation. 



^If we do not take the classical saddle point, both boundary conditions will matter. 
^For an analogous explicit derivation, see Appendix A of Ref. ES] 
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To consider fluctuations, following Litini and Manuel, we shall now define 

kA = if) + Sf (51) 

and 

A= {A) + a (52) 

where (/) and (A) denote the functional average of f^^A and A in the path 
integral for the effective action (see Eq. (|iO|)). Note that by definition {6f) 
and (a) both vanish. Furthermore, one has 

F'"' = F^"" + f"" , (53) 

where 

F"'' = d>'{AY^''-d%AY^'' + gr^^AY^^{AY^^ 
fi^u,a ^ {D^'ay-{D''a^y + gr^''a^''^a'''^ (54) 

with the covariant derivative D^ defined to be d^ — ig (A^). 
Substituting Eqs. (|51| ) and ( ^21 ) into Eqs. (|49| ) we obtain 

P^ h,{f^A^,p,Q)) -2gTT{QFn-^{kAi^,P,Q))] = {ri) + {0 ■ (55) 

In the preceding equation, (r/) and {() are given by 

iv) = 9P'Q''{f;.d;6f) 

(0 = gQT"'p'{alal)d;U)+gr"'Q'p'{aldQa5f), (56) 

and 

5m(4a(x,p,Q)) = (d,-gr'^Q^A'^dQ.) {fU^^p^Q)). (57) 

Performing the same substitution in Eq. (0), we find 

D.F^'^''^ + (JflTct) = {J'^n , (58) 

with (Jfluct) defined by the relation 
{rn = 9 r'' « {{{D^ay - {D'^dr) + g r'^ a^Xe) + {D^a^a'^f^) . (59) 

Proceeding from the preceding set of equations, Litim and Manuel have 
performed an extensive analysis, in weak coupling, of these equations and 
their solutions. For instance, neglecting the fluctuation terms introduced 



above, and expanding single particle distribution about the equilibrium Bose 
or Fermi distribution (/cquii), 

{fix,p,Q)) = {f!:Zii{Po))+9{f^'\x,p,Q)) + ■■■, (60) 

one can reconstruct the hard thermal loop effective action, as previously 
shown by Kelly et al. Keeping the ffuctuation terms, and similarly expanding 
(77), {() and (Jfluct); one obtains the Boltzmann-Langevin equation contain- 



ing both noise and collision terms that was first derived by Bodeker |]T7 . 

We will not perform such an analysis in this section, and wish to direct 
interested readers to the papers mentioned. The objective of this section was 
to demonstrate that the starting point for such analyses can be obtained, after 
clear approximations, from the many body world line formalism developed in 
this paper. It is hoped that this formalism will allow us to understand better 
the validity and the limitations of the various approximations employed, and 
how one may go beyond these approximations. 

In particular, it will be useful to better understand the following: a) 
when does the saddle point approximation for the quantum mechanical path 
integral break down? Presumably, in this case, the quantum problem can 
be reformulated in terms of Wigner functions, b) It appears that in this 
formalism one may be able to distinguish between thermal and quantum 
fluctuations. How is this achieved in practice? c) The collision-less Vlasov 
equations can be exponentiated to obtain the hard thermal loop effective ac- 
tion. The addition of damping and noise terms, as Bodeker has emphasized, 
implies that we are going beyond hard thermal loops. Does consistency then 
require that one keep the sub-leading terms in Eq. (|5^)? d) What is the 
form of the functional G[^]? Is a renormalization group treatment, wherein 
G encapsulates information about non-trivial correlations feasible? We hope 
to address some of these issues, in the finite temperature context, in a future 



work 1^1 



6 Discussion 

The world line formalism in quantum field theory, going as far back as the 
early works of Feynman and Schwinger, provides a direct representation of 
quantum field theory in terms of world lines of point particles. Recently, it 
has become popular as a powerful method to compute Feynman diagrams 
in theories with internal symmetries such as spin and color. These internal 
symmetries are treated using classical Grassmanian variables. 

In the many body problem, the world line formalism leads to an intuitive 
picture consistent with our ideas about classical kinetic theory-that of phase 
space averages over the world lines of classical particles. We have shown in 

19 



this paper that classical transport theory can be obtained form the many 
body world line formalism. As discussed by several authors, this classical 
transport theory, in the Vlasov limit sums the leading class of Feynman di- 
agrams for external momenta p ~ gT- called "hard thermal loops". The 
inclusion of noise and damping terms to give the Boltzmann-Langevin equa- 
tion, goes beyond hard thermal loops and allows one to address physics at 
the softer scale p ~ g^T. The connections of the imaginary time world line 
formalism to classical transport theory is fleshed out for the specific example 
of a 0^ theory in appendix A. 

The potential advantage of the formalism developed in this paper is that 
one has a formalism at the level of the effective action. It may therefore 
facilitate computations of effects beyond those of classical transport theory. 
Its form suggests that it is amenable to a "dynamical renormalization group" 
treatment, as has been discussed previously by Boyanovsky, de Vega, and 
collaborators |^. (See also very recent work, the many body extension of 



which may lead to a picture similar to ours [^.) Similar ideas have been used 
previously in small x physics to develop a Wilsonian renormalization group 



for wee partons |3J, ^. Indeed, one can draw a formal analogy between the 
P's and Q's in classical transport theory, and the classical sources in the small 
X effective action. Despite the formal similarity, one should caution that the 
physics in the two cases is quite different. Unlike classical transport theory, 
small X physics is dominated by light cone singularities. One consequence is 
that the light cone sources are static, not dynamic like the P's and the Q's. 
It is therefore not obvious that the renormalization group ideas of small x 
physics can be transported 'lock, stock and barrel' over to the physics of hot 
QCD. 

Nevertheless, there are connections that are worth exploring further. A. 
H. Mueller has shown recently that the classical effective action, when applied 



to nuclear collisions |47], leads to the Boltzmann equation [48]. Also, recently 
Elmfors, Hansson, and Zahed have shown that the hard thermal loop free 
energy could be related directly, via the relativistic virial expansion, to the 



S-matrix elements of high energy scattering EI 



We should point out that we have glossed over many subtleties in the 
world line formalism. In particular, how the closed time path of quantum 
field theory manifests itself in the first quantized formalism needs to be better 
understood. Whether one can separate statistical and quantum fluctuations, 
is another topic that should be studied further. Another issue that needs to 
be addressed is whether one can go beyond the one loop approximation in 
the world line formalism. This has been achieved for the vacuum case by 
Schubert and collaborators p9[, and one should be able to extend it in a 
similar fashion to the many body context. Finally, there are a large number 
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of physical problems that should be hopefully simpler to compute in the 
many body world line formalism. 
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A The effective action in the imaginary time 
world line formalism 

In this appendix we derive the leading term in the finite temperature effective 
action for a 0^ theory, using directly the one loop world line formalism at 
imaginary time. 

Starting from the full imaginary time action. 



S 



d X 



V + ^' 



(61) 



where / d'ho = Jq dr J d^x and (3 the inverse temperature T, we split the field 
into a "soft" and "hard" part, = 0^ + 0/1. In the one loop approximation, 
the action Eq. ( [HT| ) becomes. 



Si- 



loop 



d X 



-]-^sd^s + Mt + \<Ph {-d^ + 12A0; 



(62) 



Integrating out the hard part, we obtain the effective action for the soft field, 
Spff = I d'^x 



--0,920, + A0 



ii 

2 



+ ^LogDet i-d' + 12A0^) = 5^^) + SJS 



(63) 
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To evaluate the remaining determinant, we use the world line approach 



at finite temperature (and imaginary time), developed in |49], and write the 
second term in Eq. (Bl) as 



5S' = -\ E iih^-^''"" 



X 



Pt 



T>x exp 



|^*rfti|^ + 12A0^(x(ti)) 



The boundary condition Pt in the path integral is given by 

Xf,{t) = x^(0) + nP6f,4. , 



(64) 



(65) 



and for simplicity we have put the einbein e = 2. 

The term corresponding to the blackbody radiation is easily obtained by 
taking A ^ and making the substitution 



where now u^{t) = u^{0) = and 2;^ is ti-independent. We obtain. 



(66) 



si^{\ = o) 



\jd'^t Tr""^"^'"/^^"^""" 



Vae-^>'^^'''/') 



327r2 



d'z E 



dt 
— ( 
t 



-n2/32/4t 



(67) 



where we have used the explicit result for a gaussian path integral. The finite 
temperature part of Eq. (^) is given by 

1 p ^^^^ POO 



rpii „ 00 1 

7r2 j ^^ n^ 



2rpi 



Tl'T 

"90" 



ds 



(68) 



which is just the result for a noninteracting gas of scalar particles. The 
corresponding result for a free gluon gas is a straightforward generalization 
of the above procedure, using Eq. (^ to represent the internal degrees of 
freedom. 

To obtain the first nontrivial result, we calculate the corresponding cur- 
rent. 



j{y) 



6S, 



(2) 
cff 



S<Ps{y) 



i2\My) E Ht f'^^~^^'''"U ^ 



X 
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exp 



dti 



x^ 



+ 12A(/)2(x(ti)) 



dt25\y - x{t2)) 



(69) 



where the delta function is periodic in the imaginary time, with period j3. 
Note that this expression is nothing but the average value of the single parti- 
cle current, in the background field 0^. Thus, it corresponds to the imaginary 
time formulation of Eq. (^), although for a 0^ theory in the present case. 

To lowest order, we may neglect the A-dependence in the world line action. 
Using also the substitution in Eq. (|5BD, the current becomes, 

~^ <-^ dt 



j{y) 



J 



I 



n=-oo-^0 t J 

me-IoM^'/^) I d^z (f dt26\y - uih) ~z- n^t^jt) 



Now, since 



d^z (J dt26\y - u{t2) -z- npt2/t) 
the current in Eq. (|70|) becomes 

'■(^^^-^ - —My) t I 



j^'\y) 



°° dt 



dU 



-n2/32/4t 



47r2 



t' 



and the finite temperature part is 

6AT2 



j^'Ky) 



TX^ 



CO 1 

^hi) E - = ^^Vs(i/) . 



n=\ 



n^ 



(70) 
(71) 

(72) 
(73) 



lYiom. the first part of Eq. (|63D, the blackbody radiation in Eq. (68) and 
the current in Eq. (|73D , we finally get the lowest order form of the effective 
action for the soft field. 



5'. 



eflF 



d'^Z 



-dsd (f)s + A0; 



90 



(74) 



The leading term in the effective action corresponds to nothing but the 
thermal mass for the soft field, as depicted in Fig. la. In principle the thermal 
mass also depends on the infrared cutoff, i.e. the cutoff that characterizes the 
separation between the hard and soft fields, as discussed in |^. However, 
this term is subleading compared to T^ and has been neglected in the present 
treatment. 

Within the one loop approximation, higher powers of the soft field come 
from diagrams of the form shown in Fig. lb, and in general there are of 
course also corrections from higher loop diagrams. 
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soft field 
hard field 



a: b: 

Fig. 1. a-.Thermal mass for the soft field. 

b: Higher corrections to the effective action 

Let us now briefiy compare the above world line formulation with the 
classical kinetic theory, as studied in |5^, |5l|]. Within the kinetic theory we 



would write the current as 

jkin(l/) = 12A0(i/)y'rfP/(i/,p), (75) 

where the integration measure dP enforces positivity of the energy and an 
on-shell constraint, 

dP = ^^2e{po)6{pl-p'), (76) 

and f{y,p) is a solution to the coUisionless Boltzmann equation, 

(^|;^t|:)^(^-^)-(^-^'"+'^-«")/fe.rt-o, (77) 

with r the proper time. Comparing Eq. ([75| ) with the world line formulation, 
Eq. (|5^, we then have to identify 

JdPf{y,p) = { J dr5\y - xir))) = 2f:J^^ J J Vpe-I>'^^'x 



nee' ■I'o '^<l(*V4+12A</,2(x(ti))) 
Pt 



dt25\y - x{t2)) 





(78) 

where Pt is the boundary condition given in Eq. (pSf). 

In the weak coupling limit, f{y,p) can be expanded in a power series in 
the coupling constant. A, 

f{y,p) = My,p) + \h{y.v) + O(A^) . (79) 

To lowest order, the kinetic theory then gives 

jl^liy) = 12X<Piy)JdPUy,p). (80) 
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With /o equal to the equihbrium Bose-Einstein distribution ne = (e^^° — 
1)~^, this reproduces the lowest order correction from the one-loop effective 
action, Eq. (^), as expected. 

We can also compare the next correction (Fig. lb with four soft fields), 
in the limit of zero external four-momentum. Collecting the same powers of 
A, by using p^ ~ ^'^ij.'PI (which also follows from the world line action), the 



Boltzmann equation gives the solution for fi{y,p) as |^, ^, ^ 



fdPMy,p)r^<P'iy) f dP-f^ r^ ^'iy)j , (81) 

J J po apo A 

and the corresponding term for the current becomes, 

jS(z/)«AV(z/)^, (82) 

where A is the infrared cutoff on the momentum of the hard particles. 

In the world line formulation, instead of a direct infrared cutoff we will 
assume a mass m for the scalar field. Neglecting any variation of the external 
soft field, the one-loop effective action gives a contribution 

(83) 

with KQ{m(3n) the modified Bessel function of the second kind. For m(3 — > 0, 
the leading term comes from 

Y: MmPn) = — + 0(log[m/3]) , (84) 

n=l ^"^ 

which gives a current in agreement with Eq. (|82|) , provided we take m ~ A. 
This agreement between the two different methods indicates that the con- 
jectured identity in Eq. (|78D indeed is true, and gives further evidence for 
the close relationship between the classical kinetic theory and the one-loop 
effective action. In particular, Eq. (|78| ) demonstrates how the thermal aver- 
age of the microscopic one-particle distribution, J dT6^{y — x{t)), naturally 
appears in the one-loop effective action, when written as a point particle 
path integral. It also incorporates all the contributions from the soft field, 
which can then be expanded in a power series in A, just like f{y,p). Thus, 
it provides an explanation for the apparent success of using a classical trans- 
port theory for the hard particles to obtain the effective action for the soft 
modes. 
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